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Abstract 

The paramagnetic insulator to ferromagnetic metal transition in Lan- 
thanum manganites and the associated magnetoresistive phenomena is 
treated by considering the localization due to random hopping induced by 
slowly fluctuating spin configurations and electron-electron interactions. 
The transition temperature and its variation with composition is derived. 
The primary effect of the magnetic field on transport is to alter the local- 
ization length, an effect which is enhanced as the magnetic susceptibility 
increases. Expressions for the conductivity, its variation with magnetic 
field, and its connection with magnetic susceptibility in the paramagnetic 
phase are given and can be tested with further experiments. 
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I. INTRODUCTION 



Interest in mixed-valent compounds of the form A^txB^+Mn^+Mn^Oa, whose fasci- 
nating properties were discovered 1 about fifty year ago, has revived recently. The phase 
diagram as a function of x for various trivalent atoms A and divalent atoms B and the 
magneto-transport properties has been determined. 2 ' 3 For 0.1 < x < 0.3 an insulating or 
very high resistance phase at high temperatures is followed at lower temperatures with 
a metallic phase accompanied by ferromagnetic order. The basic physical point, 4 made 
by Zener in 1951, is that the elementary electron conduction process, given that valence 
states other than 3 + and 4 + are disallowed due to strong local correlations, in which two 
Mn ions of different valence interchange their valence states: Mn 3+ Mn 4+ — > Mn 4+ Mn 3+ , 
is proportional to the one-electron transfer integral t only when the initial and final states 
are degenerate. This requires that the Hund's rule coupled spins of Mn 3+ (S = 3/2) and 
of Mn 4+ (S = 2) are aligned. Otherwise, the transfer integral is of 0(t 2 /J) where J is 
the Hund's rule coupling energy, which for Mn 3+ is several eV and assumed much larger 
than t. An appropriate model Hamiltonian is 

<U> i i 

where the electrons hop only between nearest neighbor sites of different valences. For 
| J | t the spin of the 'conduction' electron is always parallel to the local Si. Projecting 
to such states 

H cff = ^t(S i ,(S j + < r) max )c 1 + c j , (2) 

ij 

i.e. the transfer integral is a function of the relative orientation of the Mn 4+ (S; = 3/2) 
spin and Mn 3+ (S = 2), i.e. (Sj + o") ma x- 

In a pair- wise hopping process the conserved value of the spin is | S; + Sj + cr | , so we 
expect the hopping to depend on this function. Semiclassically the angle 9 between two 
spins Si and S 2 is given by 

9 |Si + S 2 | , 
C0S 2 = ^S— ' (3) 

so that semiclassically the effective hopping Hamiltonian is 



H e fT = cos ^ c t c y 



(4) 



This result was first derived by Anderson and Hasegawa. 5 If one considers pairwise hop- 
ping only, one can specify the relative orientation by 9-^ alone. More generally, 6 one must 
introduce also the azimuthal angle 0ij to specify the relative orientations of spins at i and 
j so that a more appropriate form is 



0ij is a Berry phase. This may have some interesting consequences which we hope to 
discuss in the future. 

Since ferromagnetism is accompanied in these materials by metallicity, one may expect 
that magnetic polarization by an external field will cause a large drop in resistivity. The 
discovery 2,3 of such a large magneto-resistance has led to much interest in these materials. 

In the next section, I describe how an adiabatic approximation may be introduced 
by first freezing the spin configuration and calculating the electronic states. It is argued 
that the chemical potential passes from a region of localized states to delocalized states 
as x is increased. Spin polaron effects as correction to the adiabatic approximation are 
considered next as are electron-electron interaction effects which are expected to produce 
a gap in the excitation spectra and in one-electron states near the chemical potential. In 
Section III I calculate the variation of the ferromagnetic transition temperature with x. 



I am primarily interested in the paramagnetic regime where the spins {Si} are ran- 
domly oriented and fluctuate at frequencies related only to the temperature keT. Any 
spin-spin correlation only reduces the characteristic frequency. Further assume that 
k B T <C t. Then the core-spin fluctuations may be treated in the adiabatic approximation, 
in which we first freeze the spins in a random configuration {Si} and calculate the elec- 
tronic states -?/>n({Si}) in the random configuration. We can then perturbatively couple the 
thermal fluctuations of spins to find corrections to the electronic wave functions. Simulta- 
neously the potential energy of the spin configurations is given by ^ n E n ({Sj})f(E n ) where 



H eff = ^te^cosf c+ Cj . 



(5) 



II. GENERAL CONSIDERATIONS 
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f is the Fermi function. This procedure is quite analogous to the Born-Oppenheimer ap- 
proximation in electron-phonon problems, except that no independent inertia exists for 
the spins. One can also investigate the paramagnetic to ferromagnetic transition by this 
method. 

Consider then first a random configuration of spins. The transfer integrals for con- 
duction electrons are then random variables varying from to |t|. The distribution of 
the transfer integral will be derived in Section III below. This randomness was treated 
recently 7 in a dynamic generalization 8 of the coherent potential approximation (CPA) 
to calculate the equilibrium and transport properties. In another paper 9 it was treated 
perturbatively to calculate the resistivity with the conclusion that the insulating behav- 
ior in the paramagnetic phase is inexplicable in purely electronic models and suggesting 
electron-phonon interactions. This line of thinking has led to the suggestion 10 that a 
Jahn- Teller distortion is responsible for the behavior. The phenomena observed in the 
manganites occurs also in mixed-valent TmSe x Tei_ x compounds. 11 No question of Jahn- 
Teller distortion arises there. 

LaMn0 3 is an insulator in which the Mn — 6 octahedra is Jahn- Teller distorted. I 
think it is not correct to conclude that the insulating behavior is due to the Jahn- Teller 
distortion. Most transition metal perovskite compounds are insulators, even LaTi03 
where effective electron-electron interactions on Ti are smaller than in Mn due to more 
effective S-electron screening. These are all Mott insulators. Note that the derivation 
of the double exchange interaction assumes correctly that the direct electron-electron 
interaction is even larger than the exchange interaction J. If symmetry allows a Jahn- 
Teller distortion in the insulating phase, such a distortion parasitically does occur. The 
conditions for Jahn- Teller effects in metals are much more stringent and in any case not 
required by symmetry and not expected to lead to insulating phase when the metallic 
bandwidth is 0(2eV). One expects that the Jahn- Teller distortion of LaMn0 3 will appear 
in reduced magnitudes even on addition of Sr in the insulating phase and decrease in 
the low temperature metallic phase. More important than the Jahn- Teller distortions for 
quantitative purposes are the breathing mode polaronic distortions which are different 
around Mn 3+ and Mn 4+ because of the large difference in ionic size. 
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In a model for random hopping, one expects electronic states to be localized at least 
for some energies. CPA is a (self-consistent) single-site approximation while localization 
arises from interference among scatters. While the CPA is a very good approximation 
for many purposes, it cannot capture the physics of localization. I believe the transport 
properties in this paramagnetic regime for most of the range of composition as well as 
the effect of magnetic field on them can be understood only when the localization of the 
single particle states is considered. 

The present problem is one of off-diagonal randomness as opposed to the Anderson 
model 12 for disorder. This has not been so thoroughly investigated. One approximate 
but reliable calculation of the nature of the wave functions in models of off-diagonal 
randomness is due to Economu and Antoniu. 13 Their results are sketched in fig. (1). 
For a model in which the hopping matrix elements ty are randomly distributed over a 
semi-circular distribution with mean t and edges at t =F ti they find that for ti ~ t , 
which is relevant for our case, states in the energy region ||W| < E < |W| are localized 
while these in the middle are extended. W is the bandwidth of the bulk which CPA gives 
correctly as about 0.7 of the total bandwidth. 

The localization length in a problem of off-diagonal disorder is expected to diverge 
both at the band-edges and at the mobility edges. The latter is familiar, and as in 
models of diagonal disorder. The former arises because states at the band-edges are the 
exponentially rare states which travel through the crystal through routes with identical 
hopping. 

Unlike Ref. (13) the distribution of disorder in the present problem is not symmetrical. 
It follows from eq. (14) below that the distribution P(t) increases linearly with t with 
a cutoff at the maximum value. The higher moments of the distribution are therefore 
larger than in Ref. (13). So the region of localized states is expected to be larger than 
in Fig. (1). Moreover strong correlations in the present problem (in the presence of 
disorder) further favor localization. The effective fluctuations in t are further increased 
by different polaronic renormalizations of t and a Mn 3+ ion and an Mn 4+ ion. 

We may safely assume that in the frozen spin approximation, states at the chemical 
potential are localized for low dopings and that a mobility edge occurs so that for higher 
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dopings the states at the chemical potential are delocalized although strongly scattering. 
CPA should be a good approximation to calculate properties in the second regime. 



Spin-polaron Formation: First Correction to the Adiabatic Approximation 

The spins do not have any inertia. The leading correction to the adiabatic approx- 
imation then amounts to (i) minimizing the free-energy F({Sj}) with respect to spin 
configurations {S;} and (ii) calculating scattering between electronic states due to fluctu- 
ations about the new configuration. The first part is the same as considering spin-polaron 
formation. 

We can easily make the first correction to the frozen random spin-configuration by 
the formation of spin polaron around an electronic state which otherwise would be highly 
localized. Let there be P lattice sites which are spin-polarized so that an electron can 
hop freely between them to lower its kinetic energy. This is opposed by the entropy lost 
by the spins. So P is determined by minimizing 

E(P) = -W + ^ + PkT In (2S + 1) (6) 

where a is a numerical factor of 0(1), so that 

P » (at /T) 3 / 5 (7) 

This is an elementary generalization to finite temperature of Nagaoka's result 14 for such 
models that, at T = 0, a single carrier will induce ferromagnetic order. In the manganites 
the effective t (renormalized for effect other than spin fluctuations, phonons for instance) 
is about 0.2 eV. So at room temperatures the size of the spin-polaron is only a few bonds. 
So the spin-polaron effect appears not to invalidate the adiabatic approximations over 
the bulk of the band in the paramagnetic range of temperatures of interest. 

A prediction following from these considerations is that the effective moment in sus- 
ceptibility measurements even at fairly high temperatures will be larger than that of the 
appropriate average of Mn 3+ and Mn 4+ moments. At a concentration x of the spin Si 
and the rest S2, one expects 

S c 2 ff = x(Si + PS 2 )(Si + PS 2 + 1) + (1 - x - Px) S 2 (S 2 + 1) (8) 
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The estimate of spin-polaron effects would be quite altered for states in the tails of Fig. 
(1) or those near the mobility edge which have a large localization length. The polaron 
effects would have a major effect only on the surface of such states leaving the bulk of 
it relatively unaffected. So Eq. (8) is an overestimate. Moreover (8) is merely the first 
sign of fluctuations above T c which increase the susceptibility over the Curie-law. 

The adiabatic assumption made here is accompanied by the ergodic hypothesis. The 
spin-configurations change locally changing the localization length of the electronic state 
in a given small region. But averaged over the sample the distribution of the energies 
and localization lengths of the states stays the same. This process does not invalidate 
the adiabatic approximation unless the temperature is so high that several states of 
energy < 0(kT) larger than the chemical potential have spatial overlap with states at 
the chemical potential. This stability argument gets stronger due to the role of electron- 
electron interactions in depleting states near the chemical potential discussed below. 

The last paragraph is only a plausibility argument. It is a very interesting unsolved 
theoretical problem to ask for the frequency scale up to which disorder can vary without 
delocalizing states. 



The role of electron-electron interactions is always more important for localized states 
where the kinetic energy has been quenched than for Bloch states. Efros and Shklovskii 15 
have given convincing arguments that Coulomb interactions create a pseudo-gap at the 
chemical potential if it lies below the mobility edge. Therefore due to electron-electron 
interactions the density of one-particle states is modified from Fig. (I) to that depicted 
in Fig. (2). In three dimensions, the density of states in the localized regime per unit 
volume near the chemical potential is of the form 



where a is a numerical constant and k is the dielectric constant. g(e) is independent 
of the localization length near the chemical potential. Equation (9) is arrived at by 



Electron-electron Interactions 




(9) 
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considering stability of one particle excitations under excitonic renormalizations due to 
the Coulomb interactions. Corrections to it due to multiparticle excitations have not yet 
been established conclusively. 

The best way to experimentally test the localization idea and (9) is a tunneling 
measurement. Density of states of the form (9) have however already been observed in 
photoemission experiments. 16,17 The key point is that the density of states is zero at the 
chemical potential (with T/A corrections) independent of doping if the chemical potential 
is below the mobility edge and the experiment is done in the paramagnetic regime. 

Resistivity 

The conductivity of the localized states is expected to be of the variable range hopping 
form: 

a~ exp [-(To/T) 1 / 2 ] (10) 

with T ~ e 2 /K,£, for T <C T . Here I is the localization length for states near the 
chemical potential. As the effective disorder is decreased by applying a magnetic field, I 
increases leading to a decrease in the resistivity. The relationship of the effective disorder 
and the magnetization is derived in the next section. In the paramagnetic regime, the 
leading dependence of the localization length on H is 

£(H)=£ (l + x(T)H 2 /t), (11) 

where £ depends on the electron density and t is O(t). Equation (10) may be com- 
bined with Eq. (11) to suggest that magnetoresistance plotted against the magnetic 
susceptibility should fall on the same curve for different temperatures and for different 
compositions. Similar scaling behavior should be observed for other transport properties 
as well. 

Note that the leading temperature dependence of the conductivity (or magneto- 
conductivity) exhibited in Eq. (10) gives the number of carriers participating in the 
conduction. Their mobility is a weakly temperature (and field) dependent pre-factor of 
(10). This is consistent with the recent observations of Ong et al. 18 from Hall effect 
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and magnetoresistance measurements that the dominant effect of magnetic fields is to 
increase the number of carriers. 

The frequency dependent conductivity should exhibit the effects of the Efros- 
Shklovskii gap in the paramagnetic phase. A gap is indeed observed. 19 The low energy 
<j(uj) appears roughly ~ uj. This is quite remarkable and needs more detailed considera- 
tion. 

When T ~ T, the conductivity is no longer of the form (10). Now several excited 
states of energy < 0(kT) overlap the states at the chemical potential. The mobility may 
now be calculated from 

eD . . 

"=kT < 12) 

where the diffusion constant D ks £ 2 /t; £ is the typical distance between states degenerate 
to within kT which at very high temperatures is 0(a) the lattice constant; and 1/r is 
the typical frequency of spin fluctuation ~ 0(kT). So the high temperature conductivity 
approaches a constant proportional to the carrier density. 



III. ESTIMATE OF METALLIC /FERROMAGNETIC TRANSITION 

TEMPERATURE 

We can calculate the equilibrium properties in the paramagnetic phase and an esti- 
mate of the ferromagnetic transition temperature T c by a variational calculation of the 
free-energy. Let the probability distribution of the angle between neighboring spins 9 be 
given by 



P(x) = M 2 5(l - x) + (1 - M 2 )M exp 



A , 



(13) 



where x = cos 9, M is the uniform magnetization which is finite for T > T c only 
if a magnetic field H is applied, A is a variational parameter which is a function of 
temperature and Af is a normalization factor. P(cos6 l ) for A/kT — > oo aligns all the spins 
while for A/kT — > 0, the spins are completely randomly oriented provided M 2 = 0. So it 
is expected that as temperature decreases the variational calculation will lead to A(T) /kT 
changing from to very large values near the ferromagnetic transition. To calculate the 
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energy of the electrons, we evaluate the probability distribution of the nearest neighbor 
transfer integral P(t): 

P(t) = £ J\ d(cos#)P(cos£) S(t/t - cosO/2) 

= M ^(t - t ) + (1 - M 2 )AT|exp (| - l)] 

One should now calculate the electronic density of states n(E) with P(t) and calculate the 
electronic energy in terms of n(E). n(E) can not be evaluated without a lengthy numerical 
calculation. We will instead use the following approximate method: For a rectangular 
density of states with bandwidth W, with height W" 1 so that it can accommodate at 
most one electron per atom, the energy at T = for c electrons per atom is 

£ (t) = ~c(l-c). (15) 

A good approximation for the electronic energy is 

E e = f dt P(t)e(t). (16) 

Classical approximation to the entropy of spins in the large S limit using distribution 
such as P(x) give incorrect answers at high temperatures. To calculate the entropy we 
write the quantum version of P(x) by writing 



j 2 _ 4S 2 

costf - 1 = (17) 



where S is the value of the spin of an ion (Mn 3+ spin S = 2 is the appropriate value) 
and J ranges from to 2S in integer increments. Then define the quantum analog of the 
second tern of (1) 



Z = £e-w si ( Ji - 4S ^ (18) 
j 



and 



Sio„s = -^(kTlnZ) (19) 

Strictly speaking I should use (18) for the calculation of the energy also, but the final 
answer is likely to be nearly the same and the form of P(x), equation (13) leads to P(t) 
of Eq. (14) which seems physically more transparent. 
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In my calculations I have neglected the entropy of the orbital motion of the x/unit 
cell "conduction" electrons compared to that of the 1/unit cell spins. This is a small 
correction. Similarly finite temperature correction to (3) were neglected. 

Free energy then is 

F^ E c -TS spins -M-H (20) 

We may write this as 

F = t KA, T/t )M 2 + F (A 2 , T/t )] - M • H, (21) 

so that the uniform magnetic susceptibility is given by 

X - 1 = 2t a(A,T/t ) (22) 

The ferromagnetic transition temperature is given by 

a(A 2 ,T c /t ) = (23) 

The result of the minimization of the free-energy with respect to A, not surprisingly 
(since this is a mean-field calculation), is that x(T) follows a Curie law 

X(T) « (24) 

The numerical minimization of F with respect to A gives that 

k B T c ^0.1E c F oh (c) (25) 

where E^ oh (c) is the electronic cohesive energy for the ferromagnetic case 

ELW = yc(l-c). (26) 

The numerical factor of ~ 0.1 reflects (i) that the electronic cohesive energy increases by 
only about 20% in going from the complete random spin-orientation to the ferromagnetic 
configuration and (ii) that the entropic free energy of spins at a temperature T = W is 
about a factor of 2 larger than W. 

Mattheiss has calculated the conduction electron bandwidth in the local density ap- 
proximation to be 2.5 eV. Then at x ?a 0.3, the calculated T c ps 250 which, considering 
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the crudeness of the calculation, is in the right range. A more important test is the 
relation T c ~ x(l — x) which is compared with experimental results in Fig. (2). 

CPA is a perfectly respectable approximation to calculate the energetics. Therefore 
Furukawa's calculations 7 for T c are probably better than the estimates here although the 
considerations here may be more transparent. 

IV. CONCLUDING REMARKS 

An attempt has been made to understand the equilibrium and transport properties of 
mixed-valent magnetic compounds on the basis of electronic localization due to magnetic 
disorder and the alteration of this disorder by a magnetic field. A paramagnetic insulator 
to Ferromagnetic transition in the observed range of temperatures and with the observed 
composition dependence is derived. The primary effect of a magnetic field on transport 
is to alter the localization length, an effect which is enhanced as the magnetic suscepti- 
bility increases. Scaling of the transport properties with the magnetic susceptibility and 
tunneling experiments to observe the Efros-Shkloveskii correlation induced pseduo-gap 
in the paramagnetic phase are suggested. 

In this paper the ordinary antiferromagnetic exchange between localized carriers and 
the effects due to A-B disorder have been ignored. At low temperatures their effects have 
been adequately treated by deGennes. 20 

Note Added 

Based on a preprint of this note, the authors of Ref. (10) have in yet another paper on 
the Jahn- Teller distortion idea (A. J. Millis, R. Mueller, B. I. Shraiman, preprint) pointed 
out in effect that the fluctuations in hopping estimated here using the results of Economu 
and Antoniu are small by about a factor of 2 to give the mobility edge at x pa 0.3. This 
is picayune. The point of this work is to show an important effect which systematically 
explains observations as x is varied. Estimates of the fluctuation in hopping based on 
"one-electron approximation" can easily be in error by a factor of two or more. As pointed 
out here electron correlation effects increase the tendency to localization and increase the 
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effective disorder parameter. So does the different iattice distortion (breathing mode) 
around Mn 3+ and Mn 4+ which renormalizes t. The uncertanity in disorder fluctuations 
by a factor of 2 should be compared with a general lack of knowlege of coupling constrants 
in the Jahn Teller scheme. In any case recent neutron scattering results strongly disfavor 
Jahn Teller effects as responsible for the insulating to metallic crossover (M. Marezio, 
private communications). 
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FIGURE CAPTIONS 



Fig. la Character of states as a function of their energy in a tight binding model with 
disorder in which t is the mean and ti the rms width of the transfer integral. This 
sketch is based on Economu and Antoniu, Ref. (13). 

Fig. lb Sketch of the density of states. 

Fig. 2 Density of states with Efros-Shklovskii gap. 

Fig. 3 Magnetic transition temperature vs x in La7_ x Sr x Mn03 taken from Ref. (3), com- 
pared to the result T M ~ x(l — x) from Eq. (25). 
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